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In this study, we investigated the tritrophic interactions among a persistently transmitted plant virus, 
Tomato yellow leaf curl virus (TYLCV), its insect vector, the sweetpotato whitefly Bemisia tabaci, and a 
parasitoid, Encarsia formosa Gahan, one of the most extensively used biological control agents. As an 
emerging invasive pest worldwide, the two most damaging whiteflies are B. tabaci B and Q cryptic species. 
On healthy tomato plants, parasitoid-induced mortality was significantly higher in B. tabaci B than in Q. In 
contrast, similar mortality levels of B and Q were observed on TYLCV-infected plants. A higher rate of 
parasitism was consistently observed in B, independent of the TYLCV infection. Similarly, the life history 
traits of E. formosa were influenced by both TYLCV and the two cryptic species of B. tabaci. Specifically, 
E. formosa parasitizing B had a greater adult longevity and shorter developmental time on healthy plants, 
whereas the parasitoids developing from Q has a greater adult longevity on TYLCV-infected plants. The 
emergence rate of E. formosa was unaffected by either B. tabaci cryptic species or the virus. These results 
suggest that the vector-borne pathogen can manipulate the host suitability of a parasitoid and hence the 
parasitoid-host interactions. 



The sweetpotato whitefly, Bemisia tabaci (Gennadius) (Homoptera: Aleyrodidae), attacks many economically 
important crops, including vegetables, cotton, and ornamentals. Recent studies suggest that B. tabaci is a 
species complex containing at least 24 morphologically indistinguishable cryptic species' The two most 
invasive and widespread groups of B. tabaci are the Middle East- Asia Minor I species (referred to as B hereafter) 
and the Mediterranean species (referred to as Q hereafter)' ^ The B. tabaci B was first recorded in China in the 
early i990s, and it subsequently spread throughout much of the country, but has been recently displaced by Q in 
many areas^'^ 

Although B. tabaci can damage crops directly by feeding, it causes greater damage by transmitting over LOO 
begomoviruses that can have catastrophic impacts on plant growth and survival'. Tomato yellow leaf curl virus 
(TYLCV) is a single-stranded DNA plant virus that originated in the Middle East'. TYLCV is transmitted by B. 
tabaci, the only known vector, in a circulative, persistent manner^. TYLCV was first detected in Shanghai, China 
in 2006" and then quickly spread from the Southern to Northern parts of China'. Although both B and Q can 
vector TYLCV horizontally, Q exhibited a much higher efficiency than B; as a result, Q appears to be a more 
competent viral vector than B'". The spread of TYLCV coincided with the rapid displacement of the B. tabaci B 
by Q in China'. Because both cryptic species have developed resistance to various synthetic insecticides, especially 
in Q'^'", biological control is one of the best alternatives available for the sustainable management of B. tabaci 
populations. 

Encarsia formosa Gahan (Hymenoptera: Aphelinidae) is a solitary endoparasitoid that has been used com- 
mercially for the biocontrol of whiteflies, including the sweetpotato whitefly B. tabaci and the greenhouse whitefly 
Trialeurodes vaporariorum, worldwide'*. As one of the most extensively used natural enemies, E. formosa para- 
sitizes at least fifteen species of whiteflies in eight genera'*. Specially, adult females parasitize and kill whitefly 
nymph, although the developmental time of E. formosa depends on the developmental stage of B. tabaci nymph'^ 

Bemisia tabaci Q has now become the primary whitefly cryptic species in China*-', and as noted, the rapid 
spread of Q has been accompanied by the recent outbreaks of TYLCV. Information regarding the tritrophic 
interactions among TYLCV, a persistently transmitted plant virus, B. tabaci, the only known insect vector, and E. 
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Figure 1 | Biocontrol properties of E. formosa on the sweetpotato whitefly B. tabaci. The 3"^ instar nymphs of B. tahaci B and Q cryptic species were 
reared on healthy (NV) and TYLCV-infected tomato plants (V), respectively, before submitting them to a hymenopterous parasitoid, E. formosa. 
The percent mortality (A) and the parasitism rate (B) of B. tahaci were documented. Values (means ± SE) are from three biological replications. Within 
each panel, different letters indicate significant differences among treatments (P < 0.05). 



formosa, the primary biocontrol agent for whiteflies, are grossly lack- 
ing. In this study, we are specifically investigating the impacts of this 
virus and its insect vector on the performance and life history traits of 
E. formosa. A better understanding of such interactions will, in turn, 
offer valuable insight into the long-term sustainable control of white- 
flies, with or without the occurrence of epidemic plant virus. 

Results 

Biocontrol properties of E. formosa on the sweetpotato whitefly B. 
tabaci. The mortality rate of E. formosa was significantly affected by 
the B. tahaci cryptic species (F3, 236 = 34.759, P < 0.0001), virus 
infection status (F3_ 236 = 54.068, P < 0.0001), and the interaction 
between the two (F3, 236 = 23.432, P < 0.0001). The mortality rate 
was significantly lower in Q reared on healthy plants than in B reared 



on healthy plants or in either B or Q reared on TYLCV-infected 
plants (Figure lA, Figure SI). 

The parasitism rate was significantly affected by the B. tabaci 
cryptic species (F3^ 236 = 190.939, P < 0.0001), virus infection status 
(F3, 236 = 15.896, P < 0.0001), and the interaction between the two 
(F3, 236 = 27.190, P < 0.0001). The parasitism rate was significantly 
lower in Q reared on healthy plants than in B reared on healthy plants 
or in either B or Q reared on TYLCV-infected plants. In addition, the 
parasitism rate was lower in Q reared on TYLCV-infected plants 
than in B reared on healthy or TYLCV-infected plants (Figure IB, 
Figure SI). 

The impacts of TYLCV and B. tahaci cryptic species on the Ufe 
history traits of E. formosa. The developmental time of E. formosa 
was significantly affected by the virus infection status (F3_ 236 = 
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10.960, P = 0.001) and the interaction between B. tabaci cryptic 
species and virus infection status (F3_ 236 = 11.954, P = 0.001), but 
was not significantly affected by B. tabaci cryptic species (F3 236 = 
3.640, P = 0.058). Developmental time was shorter in B reared on 
healthy plants than on TYLCV- infected plants or Q reared on either 
healthy plants or TYLCV-infected plants. But it was not significantly 
different among B reared on TYLCV-infected plants or Q reared on 
either healthy plants or TYLCV-infected plants (Figure 2A, Figure 
S2). 

The emergence rate of E. formosa, i.e., the percentage of pupae 
from which adult wasps emerged, was not significantly affected by B. 
tabaci cryptic species (F3 236 ~ 0.723, P = 0.396), virus infection 
status (F3_ 236 = 0.2 3 8, P = 0.626), or the interaction between the two 
(F3_ 236 = 2.239, P = 0.136) (Figure 2B, Figure S2). 

The longevity of E. formosa was significantly affected by the B. 
tabaci cryptic species (F3_ 236 = 6.086, P = 0.014), virus infection 
status (F3. 236 = 80.644, P < 0.0001), and the interaction between the 
two (F3, 236 = 294.839, P < 0.0001). Longevity was significantly 
shorter in Q reared on healthy plants than on TYLCV-infected plants 
or B reared on either healthy plants or TYLCV-infected plants. In 
addition, the longevity of E. formosa parasitizing B. tabaci Q from the 
TYLCV-infected tomato was significantly higher than those wasps 
parasitized whiteflies from the TYLCV-free tomato. In case of B. 
tabaci B, however, the trend was reversed (Figure 2C, Figure S2). 

Discussion 

The influence of insect vector on the performance and fitness of E. 
formosa. Although parasitism of B. tabaci by E. formosa has been 
studied previously""*, none of them has been focused on the 
potential impact of cryptic species. As the most invasive and 
destructive B. tabaci cryptic species, B and Q vary in their mating 
and feeding behavior and prefer different host plants'" "'^". On 
healthy tomato plants, E. formosa performed better (higher mor- 
tality and parasitism rates) on B. tabaci B than on Q (Figure 3). 
This is consistent with the results from a congeneric parasitoid, E. 
sophia''\ In regard to fitness, the developmental times of E. formosa 
on the 3"' instar nymphs of B. tabaci were 14.16 and 14.52 days, 
respectively for B and Q. A similar trend was documented in E. 
Sophia, in which the developmental time of the parasitoid on B. 
tabaci B and Q were not significantly different"^'. Similarly, the 
emergence rate of E. formosa was not significantly affected by the 
two cryptic species. As for adult longevity, however, E. formosa lived 
longer in B than in Q. The combined results from parasitoid's 
performance and fitness indicate that B. tabaci B is a more suitable 
host for the Encarsia spp. than Q, although the underlying mecha- 
nisms are largely unknown. The E. formosa colony was obtained 
initially from the United Kingdom and was introduced into the 
vegetable-growing areas in China to control whiteflies in the 1980s. 
E. formosa has been acclimated into the field conditions in China for 
over 30 years, and experienced the invasions of _B. tabaci B and Q in 
the 1990s and 2000s, respectively^-. By 2011, Q had displaced B and 
become the dominant cryptic B. tabaci species in most areas in 
China^'^^ The present study suggests that parasitic E. formosa is 
more effective against B. tabaci B than against Q, and this diffe- 
rence in their suitability for the parasitoid could partially contri- 
bute to the competitive displacement of B. tabaci B by Q in China. 

The impact of virus infection of the host plant on the performance 
and fitness of E. formosa. Virus induced plant responses can influ- 
ence the behavior and physiology of insect vectors, and, sometimes, 
lead to behavioral changes in the vectors that favor viral transmis- 
sion^'''^^. Although tritrophic interactions among vector-borne 
pathogens, its vector, and the parasitoid (£. formosa) add another 
layer of complexity, it is germane to speculate that vector-borne 
pathogen can exert its impact on a parasitoid through the vector. 
Multiple lines of evidence have shown that TYLCV can manipulate 



the behavior (e.g. host preference) of its exclusive vector, B. tabaci 
species complex, for the benefit of viral transmission''"'^'''^'. TYLCV- 
infected whiteflies fed more readily than uninfected individuals and 
spent more time salivating, which is essential for the spread of a 
persistently transmitted plant virus'". Field surveys indicate an 
apparent mutualistic relationship between the outbreaks of 
TYLCV and the introduction and rapid range expansion of B. 
tabaci Q\ Plant defensive compounds, specifically terpenoids, play 
a key role in mediating vector-pathogen mutualistic relationships. 
Gas chromatography-mass spectrometry (GCMS) analysis showed 
that healthy plants released significantly more terpenoids than 
TYLCV-infected plants^"". The higher level of terpenoids, a suite of 
defensive compounds from plants, however, compromised the 
fitness level of B. tabaci. Suppression of terpenoid synthesis 
pathway via gene silencing, in return, enhanced the fitness level of 
B. tabaci'". That means TYLCV can influence host-derived cues in 
ways that have positive effects on vector attraction'"'^*". Among the 
two cryptic species, B. tabaci Q preferentially settled and oviposited 
on TYLCV-infected plants rather than on healthy plants, and had a 
better overall fitness, including a greater fecundity, survival rate, and 
female body size, and a shorter developmental time''"^**. 

In this study, a differential response was observed toward virus 
infection among the two B. tabaci cryptic species. Specifically, for 
B, the virus has no apparent effects on either the performance or 
the fitness of E. formosa. Increased mortality and parasitism rates, 
as well as enhanced life history traits, however, were documented 
in E. formosa parasitizing Q from TYLCV-infected tomatoes. 
Improved parasitism or predation by the virus -infected host is 
not without precedent. Belliure et al. noticed an increase in the 
developmental time of thrips on infected host plants, which led a 
reduced risk of predation by predatory mites"''. Infected with 
Tomato spotted wilt virus, thrips larvae grew too big to be preyed 
on by the predatory mites. However, when the prey size is not the 
limiting factor, plant virus could not influence the predation risk 
(a predatory beetle can consume all sizes of thrips larvae)"^'. The 
seven-spotted ladybug, Coccinella septempunctata, consumed 
more aphids when the aphids were infected with Rhopalosiphum 
padi virus^". Also, a higher rate of parasitism and the greater 
overall population suppression by the parasitoid Aphidius cole- 
mani were observed when their aphid hosts were infected with 
Cereal yellow dwarf virus^^. 

Conclusion. A persistently transmitted virus, like TYLCV, can en- 
hance the quality of host plant to attract insect vector and increase 
the viral transmission''''. The quality of the insect vector, the host for 
parasitoids, can be improved accordingly. Higher quality of the host, 
consequently, can impact the performance of parasitoids and can, 
potentially, shift their host preference (Figure 3). Previous study 
showed that TYLCV-infected plants improved the growth, 
development, and fitness of Q but decreased the fitness of B'. In 
this study, when infected, the fitness level of E. formosa decreased 
in B, whereas both fitness and performance of the parasitic wasp were 
significantly improved in Q. 

The present work showed the improved parasitism of E. formosa 
on the host of B. tabaci Q when the virus factor was added. The 
combined results suggest that the parasitoids may contribute to the 
suppression of TYLCV in the field by disproportionately influencing 
the vectors feeding on the infected plants. In addition, for the con- 
current control of B. tabaci B and Q in the field, our results are 
encouraging in the context of integrated pest management, because 
they suggest that the current biological practice with parasitoids may 
contribute indirectly to the management of TYLCV by impacting on 
more competent virus-carrying vectors, B. tabaci Q, and conse- 
quently reducing the proportion of vectors in the population that 
are infectious. 
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Figure 2 | The impact of TYLCV and B. tabaci cryptic species on the Me history traits of E. formosa. To evaluate the potential impacts of plant virus and 
insect vector on a primary whitefly parasitoid, life history traits, including the developmental time (A), emergence rate (B), and longevity (C), of E. 
formosa were recorded. The 3"^ instar nymphs of B. tabaci B and Q cryptic species were reared on healthy (NV) and TYLCV-infected tomato plants (V), 
respectively. Values (means ± SE) are from three biological replications. Within each panel, different letters indicate significant differences among 
treatments (P < 0.05). 



SCIENTIFIC REPORTS | 4 : 5926 | DOI: 1 0. 1 038/srep05926 



4 




E. formosa pupae 

Figure 3 | Patterns of parasite behavior in relation to plant pathogen and insect vectors. In this schematic, the parasitoid is an endoparasitoid, 
E. formosa, the pathogen is a plant virus, Tomato yellow leaf curl virus (infected plant is discoloured in yellow), and the vector is the sweetpotato whitefly 
B. tabaci B and Q cryptic species. Both insect vector and plant virus have a significant influence on the performance and fitness of E. formosa. Parasitic 
wasp performed better on B. tabaci B than on Q on healthy plants. Whereas, virus induced plant responses influences the behavior and physiology of the 
insect vector, and consequently exerts its impact on the parasitoids indirectly. For B. tabaci Q, TYLCV had a positive effect on both performance and 
fitness of E. formosa. In contrast, a negative impact on the life history traits of E. formosa was observed when B was fed on the infected tomato plants. Please 
note the number of B. tabaci nymphs infected with E. formosa is purely hypothetical, implicating either positive (more) or negative (less) effects. 



Methods 

Insect and parasitoid cultures. Bemisia tabaci B population was originally collected 
from an infested cabbage field in Beijing, China in 2004", and the Q population was 
collected from a poinsettia plant in Beijing, China in 2009. Their identities were 
confirmed by random amplified polymorphic DNA polymerase chain reaction 
(RAPD-PCR)^^. Both B. tabaci populations were maintained on tomato plants in a 
growth chamber at 27 ± 1' C, 70 ± 5% RH, and a 16/8 L/D photoperiod. 

Encarsia formosa culture was a gift from the Beijing Academy of Agriculture and 
Forestry Sciences. The parasitic wasp had been maintained on B. tabaci B on tobacco 
plants for more than five generations. Upon arrival in our laboratory, the original 
culture was divided into two groups; one culture was maintained on B and the other 
was maintained on Q. 

Plant materials. Tomato (L. esculentum Mill, No. 9 Zhongza) was used as the host 
plant for B. tabaci B and Q. The tomato plants were grown in the greenhouse and 
covered with 100-mesh insect-proof screen to exclude whiteflies and other pests. 

TYLCV- infected tomato plants were obtained by Agrobacterium tumefaciens- 
mediated inoculation using a cloned TYLCV genome (GenBank accession ID: 
AM282874). The plants were inoculated at the 3 to 4 true-leaf stage. TYLCV infection 
of the inoculated tomato plants was confirmed using a set of TYLCV primers C473 
and V6T'\ Bemisia tabaci B and Q populations were maintained on healthy and 
TYLCV -infected tomato plants, respectively, when host plants reached 6 to 7 true-leaf 
stage. 

Effects of B. tabaci cryptic species and plant virus on the parasitism of E. formosa. 

The parasitoid -induced B. tabaci mortality, the percentage of B. tabaci parasitized, the 
percentage of E. formosa pupae from which adult wasps emerged, E. formosa 
longevity, and E. formosa developmental time were determined using the clip-cages 
method. 3'^'^-instar nymphs were used in this study because the development of E. 
formosa is more synchronized at this life stage'^'^^. Twenty newly emerged B. tabaci 
adults (B or Q) were placed on the abaxial surface of a selected host leaf of a tomato 
plant (infected or not infected with TYLCV); the adults and leaf were contained in a 
leaf clip-cage (one clip cage per plant). After 24 h, the adults were removed, and the 
eggs were monitored and counted daily until they developed into 3'^'^-instar nymphs. 
The number of 3"^-instar nymphs was then adjusted to 30 per clip cage. One newly 
emerged E. formosa female was then introduced into each leaf clip-cage, and each of 
the four treatments (two cryptic species and two virus conditions) was represented by 
20 replicate clip-cages. The plants with the cages were maintained in growth 
chambers at 27 ± IX, 70 ± 5% RH, and a 14/10 L/D photoperiod. After 24 h, the 
adult E. formosa was removed from each clip-cage. When parasitized, each B. tabaci 
nymph was turned into one E. formosa pupa under the above mentioned 



experimental conditions. The endpoint measurements, including parasitism rate, 
mortality, and emergence rate were determined for each cohort of 30 B. tabaci 
nymphs in their respective clip-cage after 6-10 days assay period. 

A second experiment was carried out to determine the impact of B. tabaci cryptic 
species and plant virus on E. formosa longevity. This experiment used the same 
conditions and similar materials as described in the previous paragraph. One newly 
emerged E. formosa female adult was introduced into a clip-cage on a tomato leaf 
(infected or not infected with TYLCV) infested with 30 3'''^-instar nymphs of B. tabaci 
(B or Q). Each of the four treatments (two cryptic species and two virus conditions) 
was represented by 20 replicate clip-cages. Each day until the wasp died, the wasp in 
each clip-cage was transferred to a new leaf with 30 3''''-instar nymphs. The longevity 
of each E. formosa adult was recorded. Both experiments have been independently 
replicated three times. 

Data analysis. Data were first checked for normality and transformed when necessary 
to meet the assumption of normal distribution. The data for parasitism rate, B. tabaci 
mortality rate, and E. formosa emergence rate were arcsine- square root transformed 
for analyses. The generalized linear model (GLM, two-way ANOVA) was used to test 
the effects of B. tabaci cryptic species, virus infection status of the tomato plant, and 
the interaction of cryptic species and virus infection status on the following 
parameters: B. tabaci mortality rate, the percentage of B. tabaci parasitized by 
E. formosa, the emergence rate of E. formosa adults from pupae, the longevity of 
E. formosa adults, and E. formosa developmental time. SPSS (version 13.0; SPSS Inc., 
Chicago, IL, USA) was used for statistical analyses (one-way ANOVA, Tukey's test). 
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